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Abstract 

This paper studies the coherent acceleration of ions interacting with two electrostatic waves in a 

uniform magnetic field Bq. It generalizes an earlier analysis of waves propagating perpendicularly 

to Bq to include the effect of wavenumbers along Bq. The Lie transformation technique is used to 

develop a perturbation theory describing the ion motion, and results are compared with numerical 

solutions of the complete equations of motion. Coherent energization occurs when the Doppler- 

shifted wave frequencies differ by nearly an integer multiple of the ion cyclotron frequency. When 

the difference in the parallel wavenumbers of the two waves is increased the coherent energization 

Qh' of ions is limited to a small part of the phase space. The energization of ions and its dependence 

c/3 ■ on wave parameters is discussed. 

O 
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I. INTRODUCTION 



The motion of charged particles in the presence of electromagnetic waves is a rich dynam- 
ical system that has been studied for a variety of cases. Important physical applications of 
this problem occur in laboratory and space plasmas, such as for high-temperature (collision- 
less) plasma heating and current drive and the transverse energization of ions for times short 
compared to collisional times. A particular area of interest is the nonlinear heating of ions 
(as opposed to linear mechanisms such as Landau and cyclotron damping) by electrostatic 
waves propagating through a plasma in a uniform magnetic field B . 

For a single electrostatic wave propagating across B , the stochastic heating of ions by 
waves with frequency u ^> u> ci but u ^ Noj ci (where N is an integer and u ci = qB /M is the 
ion cyclotron frequency) was studied by Karney and Bers 0, 0| . It was found that ions with 
speeds across B less than the phase velocity of the wave u/k± (that is, k^pi > uj/u ci ) exhibit 
regular motion and do not gain energy. However, for wave amplitudes above a threshold 
amplitude, the ions are stochastically heated if their speeds are inside a region with a lower 
bound near u/kj_. The stochastic "webs" generated by a single perpendicularly propagating 
wave with frequency u = Nu ci also lead to stochastic ion heating ^,0, 0|. 

For a single wave propagating obliquely to B it was found that ions could also be 
stochastically heated p, l3|. It has recently been shown that single and multiple drift-Alfven 
waves with uj < uj C i can induce stochastic ion heating which may account for certain 
experimental observations |9j. For two waves propagating obliquely to B , the threshold 
wave amplitudes needed for stochastic motion can be significantly lowered There is 

still a lower bound for the stochastic region of phase space. 

For two perpendicular waves that satisfy the resonance condition u\ — u>2 = Nujd Ram 
et al. discovered numerically that coherent (as opposed to cha otic) energization can 
bring ions from low energies into the stochastic domain. Benisti et al. 2 then showed that 
this coherent energization was described by perturbation analysis using Lie transformation 
methods. The coherent energization was also shown by Ram et al. [IH to be described by a 
multiple time scale analysis, and invoked to explain the energization of hydrogen and oxygen 
ions from Earth's upper auroral ionosphere into the magnetosphere. For two non-collinear, 
perpendicularly propagating waves, the coherent energization was found to persist as long 
as the angle between them was less than 30° (lij . 

Coherent acceleration by electrostatic waves with u)\ — 0J2 = Nu C i can only occur when 
both wave frequencies are larger than uj C i. This process is most interesting for cases where 
ions with energy well below the stochastic region {k±pi <C uj/uj ci ) are accelerated into it. 
Most of the work on coherent acceleration has focused on waves with frequencies much 
higher than u ci . In magnetic fusion experiments and in the Earth's ionosphere, lower-hybrid 
waves fit this description (c^y, ~ oj pi ^> u ci , ui h = lower-hybrid frequency, u) pi = ion plasma 
frequency) . 

In this paper we study the interaction of ions with electrostatic waves ranging in fre- 
quency from lower- hybrid frequencies down to a few multiples of The analysis of Benisti 
et al. l3| is generalized to include nonzero wavenumbers along Bq. We develop a pertur- 



bation theory using the Lie transformation method and find conditions for which coherent 
acceleration persists. We also discuss the dependence of the range of energization and period 
of coherent oscillations on wave parameters. 

The Hamiltonian formulation of the problem is given in Section II. An analytic perturba- 
tion theory for the coherent motion based on the Lie transformation technique is described 
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in Section III. Section IV compares the results for the perturbation theory with numerical 
results obtained from the complete dynamical equations. The scalings of coherent ener- 
gization and the period of oscillation, for perpendicularly propagating waves, are obtained. 
Section V discusses the case of obliquely propagating waves and compares the results with 
those for two perpendicularly propagating waves. 

II. EQUATIONS OF MOTION 

The nonrelativistic equation of motion of an ion in the presence of a uniform magnetic 
field Bq = Bqz in a plasma and interacting with two electrostatic waves is 

M— = q sin(fcj ■x-u i t + a*) + qv x B (1) 

i=l 

where $j is the electrostatic potential amplitude, ki is the wavevector, Ui is the wave fre- 
quency, and «j is the phase of the i th wave. We normalize times to the inverse of the ion 
cyclotron frequency u ci , distances to the inverse of ki x , and masses to the ion mass M. We 
restrict our attention to the case where both k^s lie in the x — z plane. Let V{ = Ui/ud and 
ei = (ujBi/uci) 2 , where ujbi = (qkl x Qi/ M) 1 / 2 is the bounce frequency in the i th wave. The 
Hamiltonian for this system is 

h(x,p,t) = l(p — A) 2 + ^eiCos(fci ■ x — Uit + ofj) (2) 

i 

where A = B xy is the vector potential, and p = mv + qA — > v + xy is the (nondimensional) 
canonical momentum. 

Since h is independent of y, we can eliminate the y degree of freedom by making a 
Galilean transformation to a frame moving in the y direction with speed p y o = v v q + Xq (the 
subscript refers to a quantity's initial value). Following Ref. [l5j], the generating function 
for the canonical transformation from (y,p y ) to (y',p' y ) is F 2 = (y — p y ot)(p' y + p y0 ). Then 
y' = y — p yQ t and p' y = p y — p y0 , so that p' y0 = 0. The transformed Hamiltonian (to within a 
constant) is 

h(x,z,p x ,p' y ,Pz,t) = \ [p 2 x +P y + (x -p y0 ) 2 +P 2 ] + ^2 / e i cos{k i ■x-u i t + a { ) (3) 

i 

Since dh' /dy' = 0, p y is independent of time so that p' y = 0. This eliminates the y' degree 
of freedom from the dynamics. Replacing x by x' = x — p y o and p z by v z gives 

h'(x', z,p x ,v z ,t) = \ (p 2 x + x' 2 + v 2 ) + ^ €i cos(k ix x' + k iz z - u{t + a { ) (4) 

i 

where a, + k ix p y0 is replaced by «j. 

In a frame moving with velocity uz the Hamiltonian remains unchanged except that the 
wave frequencies are Doppler shifted: v% — > Ui — k{ Z u. Without loss of generality we assume 
that v z q = and consider the effects of v z q on oblique propagation in Section El 
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We transform (x', p x ) to action-angle coordinates (0, J) using the generating function 
F\ = \x' 2 cot0. Then I = \{y 2 + x' 2 ) = \{y 2 + v 2 ) is the perpendicular kinetic energy, and 
= arctan(x'/f x ) = arctan(— v y /v x ) is the gyrophase. The transformed Hamiltonian is 

H (0, z,I,v z ,t) = I + \v\ + 2J e» cos(k ix p sin + A; i2 2; - z/;t + a 4 ) (5) 

i 

where p = \/2J is the ion gyroradius. 



III. PERTURBATION ANALYSIS OF COHERENT MOTION 



In general, the equations of motion obtained from (J5J) cannot be solved analytically. 
Consequently, we resort to numerical solutions to provide an insight into the dynamics of 
ions in two electrostatic waves. Figure Q shows the time evolution of p for three ions having 
the same initial p , but different 0o, interacting with two waves of frequencies vx = 40.37 and 
u 2 = 39.37, and amplitudes ex = e 2 = 4. (All the numerical solutions of ordinary differential 
equations have been carried out using the Bulirsch-Stoer algorithm described in Ref. (l6|-) 
There are two distinct kinds of motion: the slow, smooth, "coherent" oscillations at lower 
p, and the irregular, "stochastic" motion at higher p. Superimposed on the coherent motion 
are small-amplitude, high-frequency fluctuations. Figure El shows the orbits for the same 
parameters as Fig. ^except that v 2 = 39.369 and the initial conditions are different. This 
demonstrates that the coherent acceleration from low to high energies occurs only when 
Vx — v 2 is an integer. 

Our interest is to provide an analytical description of the coherent dynamics without 
going into details of the stochastic region, other than to note its existence for p rs min(z/j) 
|2| • We assume that the waves are perturbing the cyclotron motion of the ions and express 



H = H + Hx (6) 

where 

H = I + and H x = cos(/c ix psin + k iz z — vfi + a»). (7) 

i 

An approximate analytical description of the ion motion in the coherent regime is obtained 
by using the Lie perturbation technique [13, with the ordering parameter e (e ~ e\ ~ e 2 ). 
We assume that Ui ^ Z but {u\ — u 2 ) = N e Z. For Vi G Z a web structure is formed 
in phase space and has been discussed elsewhere for a single wave |4| and for two waves 
propagating across Bn ~vh. For the case of a single wave the stochastic web structure also 



has a lower bound 

From the Lie perturbation analysis (Appendix EJ the Hamiltonian that describes the 
coherent ion motion to 0(e 2 ) is 



H(<P,z,I,v z ,t)=I + ±v 2 z +H 2 (8) 
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where 



H 2 = S (I,v z ) + S-(I,v g )cos(N{<l>-t) + Ak z z + ax - a 2 ) (9) 

So = 5*02: + So z (10) 
Sqx = —^/.kix 6 ! Jm,iJ m i (11) 

2p / — ' m — Hi 

S_ = SL* + S- z (13) 

5-3, = ——. s(ki x (m- N)J' ml J m _ N2 + k 2x mJ m>1 J' m _ N 2 ) (14) 

4p(m-/ii) 

(ki x mJ mj2 J m+N2 + k 2x (m + N)J m2 J m+N ^) 



Ap(m - p 2 

n li 7 / Jm,lJm-N,2 , Jm^Jm+NA \ /-. r \ 

A/c 2 = (/c l2 — k 2z ), ^ = V{ — k iz v z , and m is summed from — oo to +oo. J m i = J m (k ix p) is the 
Bessel function of the first kind and /'(£) = df/d£. The barred coordinates are related to 
the original coordinates by a near-identity transformation: (0, z, I, v z ) = (cp, z, I, v z ) + 0(e) 
(Appendix EI). For instance, the relation between I and I is: 

I ~ I — f i — cos(m0 + - z/jt + ai) (16) 
' m — pi 



The Hamiltonian H is a generalization to oblique waves of the results obtained in |12| for 
collinear perpendicularly propagating waves. In the limit k^ z — > the above reduces to the 
description in Ref. 12 . A nonzero ai — a 2 is equivalent to a shift in the initial 0o so that, 



without loss of generality, we can set ot\ = a 2 = 0. 

Our perturbation analysis assumes there are no resonances at 0(e). Such resonances 
occur if i>i is an integer, where our present analysis breaks down. 

The explicit time dependence in H can be eliminated by transforming from to ip = — t 
using the generating function F 2 = 7(0 — t). The transformed Hamiltonian is: 

H$, z, I, v z ) = \v\ + S (I, v z ) + SL(I, v z ) cos(A^ + Ak z z) (17) 

where 1 = 1 has replaced I (I and ip are canonically conjugate). Since H does not depend 
explicitly on time it is a constant of the motion. 

Using Hamilton's equations for I and v z , we find a second constant of the motion: 

" - = (18) 



dt V N 

Thus, the system is integrable and the dynamics described by H are not stochastic. Along 
an orbit, v z is a function of I and initial conditions only: 

Ak 

v z = v z0 + ^(I-I ). (19) 
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Therefore, Sq and S- are functions just of I. Since | cosx| < 1 

H <H< H+, H ± (I) = \v\ + ^(J) ± |5_(J)| (20) 

For an initial condition with a given value of H, p varies between the two points where H 
equals H + (p) or if_(p). We refer to H± as the potential barriers, since they delimit the 
allowed and forbidden regions of phase space. 

Figure H3 shows the orbits generated by the second-order Hamiltonian (fTTj) for the same 
parameters as in Fig. ^ Our perturbation analysis accurately captures the coherent motion 
of the full system except near the stochastic region p w min(i/j) where our perturbation 
theory breaks down. Below this region, p and p differ by small fluctuations that are accounted 
for, to 0(e), by the transformation (|16j). 



IV. COHERENT MOTION FOR PERPENDICULAR WAVES 

Using the Hamiltonian (|17|) we now analyze the ion motion for two perpendicularly prop- 
agating waves. Figures n and EH show the complete and coherent motion, respectively, for two 
perpendicularly propagating waves (k iz = 0). Figure |U displays H + and H_ from (J2UJ) for 
the same parameters as in Figs. Q and El and the values of H for the three initial conditions. 
The coherent analysis correctly predicts that particle 3 in Fig. ^ will not make it into the 
chaotic regime because it is reflected by the bump in if_. 

If we multiply e% and e 2 by the same factor a then H in (fTTj) is multiplied by a 2 (note 
that for perpendicularly propagating waves, \v\ is a constant and can be eliminated from 
H). Since a rescaling of the Hamiltonian is equivalent to a rescaling of time, rescaling both 
6j's does not affect the range of motion in p but rescales the period by 1/a 2 . For e\ ~ e 2 , 
this means the period scales like 1/ef . This reflects the fact that the coherent motion is 
second-order in the wave field amplitudes. It also shows that in certain physical situations, 
at sufficiently large amplitudes, the effects of collisions on the coherent energization can be 
made negligible. 

The range of coherent motion in p scales linearly with the wave frequencies. In Fig. E] 
we plot H±/H + (^ = 0) versus £ = p/v\ for two values of v\ with N = 1. Note that the 
potential barriers do not change significantly with v\. Figure El shows that, as a function 
of Ui, the average £ m i n and £, max have a small variation (£ m j n and £ maa; are the maximum 
and minimum £ attained by an ion undergoing coherent motion and occur when the ion 
reaches the barriers H±). The average is over ions with the same initial £o — 0.4 and 
different O — (0, 0.05, . . . , l)ir. Waves with higher frequencies can therefore produce co- 
herent energization to higher energies. Since the lower bound of the stochastic region is 
p w min(i/j) — > £ sa 1, ions with the same initial £o and 0o either will or will not reach the 
stochastic region regardless of the wave frequencies. For v\ near an integer £ max is about 
20% higher than when v\ is near a half-integer (this is not shown explicitly in Fig. 

The period of oscillation in I (see Fig. |3J) can be estimated from the equation of motion 

for J: 

PjTJ 

I = =■ = NS„ sin( Nip) (21) 

dip 

An orbit's turning point typically occurs when ip = nir/N, i.e., when it hits one of the 
barriers H±. Therefore, approximately, the period of oscillation r is given by r w 27r/ (iV (?/>)), 
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where () denotes the average over one period. From the asymptotic forms of So and S- for 
v\ ~ V2 ^ 1 (Appendix EI) we find that 



H^vfh a {^) (22) 
where h a depends on Vi only through £. Then 

j - « u' 2 ^ - (23) 
* ~ 81 1 81 1 i {l6) 

Thus, t ~ v\. Waves of lower frequency accelerate ions much more rapidly than those 
with higher frequency and thus may also be made less sensitive to the effects of collisions. 
Figure E| compares this scaling with the periods of two actual orbits obtained from H . 



V. COHERENT MOTION FOR OBLIQUE WAVES 

In this section we describe the motion of ions when the waves have nonzero parallel 
wavenumber ki Z . This couples the parallel dynamics to the perpendicular motion. 

For ions with initial v z q = interacting with a single oblique wave, the motion is stochastic 
when 0|: 

\/\Jn (p)\ + y/\Jno+l(p)\ > (24) 

where n Q is the greatest integer less than v. For n ^> 1, the lower bound of the stochastic 
region is 

P - n + -^-n /3 - l.lnl /S (25) 

As for a single perpendicular wave, the lower bound in p is roughly the wave frequency and 
decreases with e. The stochastic region in v z extends from v z ~ to v z ~ 2v. 

For two waves Eq. (|19|) shows that v z changes only when /\k z ^ 0. Thus, the cases 
Ak z = and Ak z ^ lead to different dynamics and are treated separately. 



A. Equal Parallel Wavenumbers: Ak z = 

Figures |H1 and El show the time evolution of p and v z for two waves propagating at an angle 
of 45° {k iz = k ix = 1) to B . As in the case of two perpendicularly propagating waves, there 
is coherent change in p. During this coherent evolution v z has small-amplitude fluctuations 
around its initial value. In the region where the motion in p becomes stochastic so does the 
motion in v z . The stochastic region in v z agrees with the estimate given above. Since v z is 
a constant during the coherent motion, the fluctuations in v z are due to the transformation 
between v z and v z . 

The main effect of equal parallel wavenumbers is to slightly decrease the range of coherent 
motion from what it is for perpendicularly propagating waves, thus inhibiting some ions from 
reaching the stochastic region. Numerical studies indicate that \S 0x /S 0z \ and \S- X /S- Z \ are 
both unity for £ < 1 but approach as £ — > 1. This raises the bump in H_ as k iz is increased. 
Consequently, more ions are reflected by H_ and the range of coherent motion in p is slightly 
lowered. This is evident from Fig. which shows H±/H + (p = 0) for k\ z = 0.1 and 1. 
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Figure HP shows the range of motion £ m j n , ^max for different k\ z . Increasing k{ Z slightly lowers 
£,max since the enhanced bump in H- reflects more ions. Generally then, significant coherent 
energization and access to the stochastic region is obtainable with oblique waves provided 
that Ak z = 0, while the normalized k z may be large. 



B. Unequal Parallel Wavenumbers: Ak z ^ 

When the parallel wavenumbers of the two waves are different, the coherent motion of 
the ions changes drastically. In this case v z undergoes coherent motion and the term in 

H (fTTj) is no longer a constant. This limits the range in p as Ak z is increased. Figures El 
and show the time evolution of p and v z for the exact orbits obtained from (J5J) with 
k lz = 0.001 and k 2z = 0. These figures illustrate the limits in p. 

Figure [HI shows the variation of H± and \v 2 z as functions of p. For p far from p = 17, 
H + — H_ = 2\S-\ <C \v 2 z so that H + « H_. Figure IT^l shows the limitation on the range of 
coherent motion in £ for Ak z ^ 0. 

The coherent motion in v z has the effect of detuning the waves from exact resonance. 
The resonance condition for an ion with v z ^ is 

R = v 10 - v 20 - (Ak z )v z G Z (26) 



z ■ 



where v w and v 2Q are the wave frequencies in the laboratory frame, and Vi = v i0 — k iz v 
H describes the ion's motion as long as R is close to an integer. For Ak z ^ 0, v z changes 
coherently. Condition (|26|) is not satisfied for all times, and the resonant interaction becomes 
less effective. The coherent change in v z thus limits itself, which keeps R close to an integer. 
Since the coherent changes in p and v z are linked via (|19)1. the coherent change in p is also 
small. 

Consider a distribution of ions with different initial v z q interacting with two waves of 
frequencies v\q and u 2 o- For Ak z = 0, all ions will be in resonance with the waves provided 
b>io — V20 G Z. For Ak z ^ 0, the resonance condition implies that only ions with certain 
v zQ , namely 

Vza ~ ^ , neZ (27) 

are initially in resonance. As v z changes coherently, they fall out of resonance. 

This situation is analogous to the case of two perpendicularly propagating waves when 
the wave frequencies do not differ by an integer 0. Following Section IV. C of 0, the 
approximate Hamiltonian, correct to second order in wave amplitudes, that describes the 
coherent motion is 

Hoff = ~I + H (28) 
where — v 2 ) = N + Av and \Av\ 1. In this case the barriers H± are given by 

H ± = -^fl + ±v 2 z + S ±\S_\ (29) 

The first term in (|28|) limits the coherent motion, and plays a similar role to \v\. Figure ITBI 
shows the range of motion in £ as a function of Av. the largest range of coherent motion 
occurs for Av slightly different from 0, which allows —(Av/N)I to partly cancel So in 
Eq. flUD- 
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As the wave frequencies are increased, the range in Ak z and Av for which there is appre- 
ciable coherent motion becomes much narrower. Let !; a (Ak z ) be either the upper or lower 
bound of coherent motion in £ for wave frequencies v\ a and v 2a = v\ a — N . The asymptotic 
forms in Appendix |B] indicate that for two different frequencies v^b and u 2b = v\b — N, 



Suppose £ a is large for k\ < Ak z < k 2 , and that v xh = 4z/ la . Then £ 6 is large only for 
fci/64 < Ak z < ^2/64. Coherent motion occurs over a smaller range of Ak z when the wave 
frequencies are larger. Similarly, 



Figures ITB1 and Hoi demonstrate the range of coherent motion versus Ak z and Av, respectively, 
and validate the scalings in (|30|) and (|31|) with wave frequency. As the wave frequencies are 
increased, Ak z and Av must be much smaller for ions to be energized to the stochastic 
region. Hence, just as in the cases of perpendicular propagation or Ak z = 0, for nonzero but 
small Ak z energization by waves with low frequencies is more advantageous than by waves 
with high frequencies. 

VI. CONCLUSIONS 

We have shown that two electrostatic waves propagating obliquely to an ambient mag- 
netic field can coherently energize ions when their Doppler-shifted frequencies differ by a 
multiple of the ion cyclotron frequency. A second-order Hamiltonian, derived using the Lie 
perturbation technique, accurately describes the coherent motion and agrees well with nu- 
merical simulations of the complete dynamical equations. The energization of ions occurs 
regardless of the angle of wave propagation provided the parallel wavenumbers of the two 
waves are approximately equal. If the parallel wavenumbers are equal, there is no coherent 
acceleration along B but considerable stochastic energization both along and across B Q . 
Moreover, the perpendicular coherent motion is quite similar to the case of perpendicularly 
propagating waves. There is a small amount of coherent acceleration along B when the 
parallel wavenumbers differ, but this causes the resonance condition to be violated. A dif- 
ference between the parallel wavenumbers is similar to the difference between [uo\ — u^)/^ 
and the nearest integer. 

There is no threshold ion energy or wave amplitude required for the coherent acceleration. 
The change in the ion gyroradius is linear in the wave frequencies and independent of wave 
amplitude. The period of coherent motion is inversely proportional to the square of the wave 
amplitudes and is proportional to the fourth power of the wave frequency u (uj ~ u>i ~ u 2 ). 
Furthermore, the deviation from resonance Auo = uj\ — uo 2 — Ncu ci for which appreciable 
coherent acceleration occurs scales like uj~ a , while the range in Ak z = k\ z — k 2z for coherent 
motion scales like u~ 3 . This implies that for lower- frequency waves coherent ion acceleration 
is faster and less sensitive to small changes in wave parameters. 

Coherent ion energization occurs for two waves with appropriately chosen frequencies. 
An experiment is being constructed that will be able to test the theoretical predictions of 




(30) 




(31) 
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this paper j2l|. Coherent acceleration could also occur for a broadband spectrum of waves 
extending over at least two uo C i in frequency. Such a situation can occur naturally in the 
Earth's ionosphere [13]. Detailed analyses of a broad spectrum of waves, and of the effects 
of weak collisions, remain to be carried out in future work. 



Acknowledgments 

The authors thank Prof. A. Brizard for helpful discussions on the Lie transforma- 
tion technique. We also appreciate enlightening comments from Dr. D. Benisti about 
his work on this problem. We thank R. Spektor for discussing his work with us and ex- 
ploring possible experimental realizations of this process. This work was supported by DOE 
Contract DE-FG02-91ER-54109, DOE/NSF Contract DE-FG02-99ER-54555, NSF Contract 
ATM-01 14462, and Princeton University Subcontract 150-6804-1. DJS was partly supported 
by an NDSEG Graduate Fellowship. 



APPENDIX A: LIE PERTURBATION METHOD FOR TWO OBLIQUE WAVES 



We develop the Lie perturbation method following Refs. [18[ and |17[ and follow the 
notation in Section 2.5 of Ref. [l8| . 

The Lie method provides a Hamiltonian H that describes just the coherent motion, and 
a change of coordinates that accounts for the incoherent fluctuations. The physical variables 
x = (q,p) are governed by the full Hamiltonian H(x), and the new coordinates x = (q,p) 
are governed by H(x). x depends on x and a parameter e which orders the perturbation via 

Ox 

— = [x, w(x,t)} s , x(e = 0) = x (Al) 

where [f,g) x = "%2i[(df/dqi)(dg/dpi) — (9f/dpi)(dg/dqi)] is the Poisson bracket. The old 
coordinates enter only as a condition for e = 0, which ensures that the transformation for 
any w is canonical and near-identity. 

The operator T relates the representation of a physical quantity / in the two coordinate 
systems by f(x) = (Tf)(x). In particular, f(x) = x gives x = Tx. T satisfies 

^-f(x) = -T[w(x,t)J(x)} x (A2) 



H is given by 



H{x) = T~ l H{x) + T- 1 I de'Tie' f™^^ (A3) 



o 



The second term is not needed for an autonomous system. 

We expand w, H, T, and H in powers of e and equate terms at each order in e. Collecting 
terms in (|A3J) at each order in e gives equations for -uv Upon carrying out the perturbation 
expansion to second order in e, we find 

D oWl =H X -H X (A4) 
D w 2 = 2(H 2 - H 2 ) - [w u Hi + H x \ (A5) 
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D Q f = dtf + lf, Ho] is the time derivative along the unperturbed trajectories. All expressions 
here are functions of the same set of coordinates. For simplicity we use x for this purpose, 
but the final expression for H governs the evolution of x. Clearly, Hq = H$. 

For T to be a near-identity operator, w must remain small. We choose Hi in the right- 
hand side of (jA4j) and (jA5|) to eliminate any terms that would violate this condition. Such 
terms are referred to as "resonant" terms. 

For the two- wave problem, H and Hi are given in (J7|), while H \ = for % > 2. Using a 
Bessel-function identity (see p. 361 of |22j), we obtain 

oo 

Hi = 22 e iJm,i cos 1p mi (A6) 

i m=—oo 

where J m ^ = J m {kixP) and ip m i = rri(p + k iz z — vfi + a«. Then from (|A4|) 

(d t + 8^ + v z d z )wi =Hi-^2 €iJ m> i cosipmi (A7) 

The unperturbed orbits are (ft = t + (fi , z = z ,v z = 0, 1 = I (in a frame where the ion's 
initial v z q = 0). Along these orbits there are no resonant terms on the right-hand side of 
(|A7|) . so we choose Hi = 0. Then 

wi = - y~] *T — sin ip mi ( A8) 

^ m + - z/j 

Since if 2 and i/i are zero, (|A5J) leads to 

(3t + <9 + «,9,)«; 2 = 2# 2 - K ffi] (A9) 

From (|XSjl 



[iwi,£Ti]= j^- — - — (-k jx mJ mA J' nJ + k ix nJ' mi J nJ ) cos T 

2p m — fii 

Jm,iJn 

(m - \i % 



+ e ^ 2 - (-k jx mJ mti J' n • - k ix nJ' mi J nij ) cosT_ (A10) 

2p m — ^ ,J 

+ \eiejk iz k jz Jm ^ n ^ (_ cos r + + cosr_)| 



where r± = ± ^ nj -. Along the unperturbed orbits, r± = {m — ± (n — + const. 
Some terms are resonant when i = j regardless of the Vi's. Other terms are resonant when 
either 2z/j, N + = {i>i + u 2 ), or N = (u t — u 2 ) is an integer. We construct H 2 to cancel these 
terms: 

H 2 = S (I,v g ) + 6S-(I,v z )cos((i/i - f 2 )(0-t) + (feu - fez)^ + «i - a 2 ) 

+ tf+SVC^z) cos((z/i + u 2 )(4> -t) + (hz + hz)z + «i + as) / All x 
+ 22 SiSii 1 ' v z) cos(2z/j(0 - t) + 2k iz z + 2«j) 
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where 5-,5+, and Si are unity when, respectively, N, N +} and 2z/j are integers and otherwise. 
Equations (fTUjl and (fT3|l give So and S 1 -, and 

S+ = ~ Yl Upfm- u ) ( klx ( m ~ N +) J 'mA J -rn+N+,2 + k 2x mJ m)1 j'_ m+N+t2 ) 
!, , Jm,lJ-m+N +! 2 

(m-fa) 2 

+ the same with subscripts 1 and 2 switched j (A12) 

Si = — S — t A m kixJm,iJ- m + 2Ult i + ( m — ^^i)kixJ mi iJ~m+2ui,i) 

y^pym fa) 



+H4 J ? <J " m+ y } (Ais) 

(771 — yUi) J 

The coherent Hamiltonian is 

H(x,t) = H Q (x) + H 2 (x,t) (A14) 
Using ip = (p — t as the coordinate conjugate to J, the transformed Hamiltonian is 

H$,z,I,v z ) = lv 2 z + H 2 (A15) 

if is a constant of the motion. 

When only one resonance condition is satisfied, we find a second constant of the motion 
besides if, which relates J and v z (when only N is an integer, this constant is v z — (Ak z /N)I). 
The dynamical system described by (jA15|) is thus completely integrable. When any two 
resonance conditions are satisfied it is easy to see that v\ and v 2 must both be half-integers. 
Then all four resonance conditions are satisfied. It does not appear that, in this case, there 
exists a second constant of the motion. The dynamics described by H could be stochastic. 

To find the transformation relating x and x, we expand (jA2|) and use x = Tx. To first 
order in e we obtain 

x = Tx^x-e[w 1 (x,t),x} x + 0(e 2 ) (A16) 

As desired, the coordinate change is near-identity. The relation between I and I is given in 
(HEJ). 

APPENDIX B: ASYMPTOTIC FORMS FOR 5 AND 5 

Here we derive the asymptotic forms for the terms in H, given in ([17)1 . using results in 
Refs. H H H 0. For k lx = k 2x = 1, let 



Sox = ^e 2 s x (p,/ii,0) (Bl) 
i=i 

S- x = exe 2 (s x (p, fa, N) + s x (p, /i 2 , -N)) (B2) 

2 

Soz = 5>^*(P>/^0) (B3) 

i=i 

S- z = e 1 e 2 k lz k 2z (s z (p,fa,N) + s z (p,fi 2 ,-N)) (B4) 
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where 

7T 

s x (p,(j,,n) = (-) n -csc/i7r(J At+ iJ_ (M+ i )+n - J M _iJ_ (At _ 1)+n ) (B5) 
o 

7T d 

s z (p,p,n) = (-) n+1 -—[cscp7iJ^ IM+n } (B6) 

and J M = Jpip). 

Bessel functions of negative order are replaced with 

J_ m ps -sin(/i7r)Y^ (B7) 

where is the Bessel function of the second kind. Using the asymptotic forms for J M and 



Yfi |22| . we obtain 

a(/i,p,n) = J M (p)Vn(p) ~ ^ (B8) 

/? = — n) tanhaotanha n ) -1 ^ 2 (B9) 

71 

7 = //(tanhao ~ ^o) ~ (/•* — n)(tanha n — a n ) (BIO) 

sech a n = — - — (Bll) 
li — n 

For N + 1 <C p all the n/ /x's are small, and to leading order in n//iwe find 

(3 ps -^(l-^ 2 )- 1 ^ ( B12 ) 

7 ps -na (B13) 

a « .J.^^) (B14) 

a (e,n) ee {l-er 112 i —^== ) (B15) 



Thus, 



7T 9 

s 2 PS -—a(p,p,n) (B16) 



4 9/i 
/(g,n) 
A* 3 



Similarly, 



(B17) 



/(£,n) = ifa + ^) (CIS) 



4 V ^ 



7T 

Sx ~ -(a(ii + l,p,n)-a(ii-l,p,n)) (B19) 

ps |( ai ( e )_ ai (_ e )) (B20) 

* l(c) " Mi + (B21) 
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where e = is small for ji ^> 1. Expanding <7i to leading order in e, 



Thus, 

A' 2 

So and S- both scale like 1/z/j 3 when z/! ~ z/ 2 - 



S;r ~ ^(o) (B22) 
= (B23) 



s x « ~ -1 (B24) 
fi 2 
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FIG. 1: p versus t for three ions interacting with two perpendicular waves from the full Hamiltonian 
H ©. Quantities in all figures are given in terms of the normalized units defined in the text. The 
initial po = 15.95 (£o = Po/^i = 0.4) for all three ions while their phases are 4>o = (— 0.3, 0.2, 0A)ir 
for ions labelled 1, 2, and 3, respectively. The parameters for the two waves are: ei = e 2 = 4, 
hx = k 2x = 1, hz = k 2z =0 3 i/i = 40.37, and z/ 2 = v\ - 1. 
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FIG. 2: p versus t for the same parameters as in Fig. ^ except that po 
U2 = Ul - 1.001. 
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FIG. 4: and iJ_ versus /9 for the same parameters as in Fig. ^ The initial values of H for the 
three ions in Fig. Qare marked by the open circles. 
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FIG. 6: Average £ m j n and £ m ax versus v\ for perpendicularly propagating waves based on the 
barriers H±. Parameters are as in Fig. |S] except that v\ = (3.37, 4.37, 80.37), £o = 0.4, and the 
average is over 4>q = (0,0.05, l)ir. 
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FIG. 7: The period of coherent oscillation r versus v\. The other wave parameters are the same 
as in Fig.^ The initial conditions are £o = 0.4 and cfto = (0.2, 0.6)-7r. The open circles and stars are 
the periods obtained from integrations of the dynamics given by H. The solid lines are proportional 
to uf with the constant of proportionality chosen to match the period at v\ = 40.37. The vertical 
axis is scaled so that v\ is a straight line. 
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FIG. 9: v z versus t for the same parameters as in Fig. |H1 
10 




FIG. 10: H±/H + (p = 0) versus p for k\ z = k2 Z = 0.1 (dots) and 1 (dashes). Other parameters 
are as in Fig. El The curves for k\ z = are very close to those for k\ z = 0.1. 




FIG. 11: Average £, m in and S, ma x versus k\ z for &2z = ki z from H±. v\ = 40.37 and the other 
parameters are as in Fig. HO 
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FIG. 15: Average Cmin and £, m ax versus Ak z from H± for k\ z = 0.1, k 2z = k\ z — Ak z , v\ = 10.37 
(dotted) or 40.37 (dashed), and the other wave parameters as in Fig. ^ Initial conditions and 
averaging are as in Fig. The abscissa for v\ = 40.37 has been rescaled by (10. 37/40. 37) 3 . If the 
scaling in (|3*U)) were exact, the dotted and dashed lines would coincide. 
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FIG. 16: Average (, m in and £, ma x versus Av from H± for v% = v\ — 1 — Ai/, v\ = 10.37 (dotted) or 
40.37 (dashed) and the other wave parameters as in Fig. ^ Initial conditions and averaging are as 
in Fig. H3 The abscissa for v\ = 40.37 has been rescaled by (10. 37/40. 37) 4 . If the scaling in (|31|) 
were exact, the dotted and dashed lines would coincide. 
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